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The efficient electrocatalysts for many heterogeneous catalytic processes in energy conversion and storage 
systems must possess necessary surface active sites. Here we identify, from X-ray photoelectron 
spectroscopy and density functional theory calculations, that controlling charge density redistribution via 
the atomic-scale incorporation of heteroatoms is paramount to import surface active sites. We engineer the 
deterministic nitrogen atoms inserting the bulk material to preferentially expose active sites to turn the 
inactive material into a sufficient electrocatalyst. The excellent electrocatalytic activity of N-In 2 0 3 
nanocrystals leads to higher performance of dye-sensitized solar cells (DSCs) than the DSCs fabricated with 
Pt. The successful strategy provides the rational design of transforming abundant materials into 
high-efficient electrocatalysts. More importantly, the exciting discovery of turning the commonly used 
transparent conductive oxide (TCO) in DSCs into counter electrode material means that except for 
decreasing the cost, the device structure and processing techniques of DSCs can be simplified in future. 



Increasing energy demands have stimulated intense research on electrocatalysts for oxygen reduction reaction, 
water reduction reaction and triiodide (I 3 ) reduction reaction etc, which are at the heart of the typical key 
renewable-energy technologies such as fuel cells 1,2 , water splitting 3,4 , dye-sensitized solar cells (DSCs) 5,6 and so 
on. Platinum (Pt) is always the suitable electrocatalyst for these reduction reactions due to its high conductivity, 
good catalytic activity and chemical stability, but its low abundance ratio and high costs precluded the large-scale 
utilizations of these renewable-energy technologies. These considerations have led to ongoing efforts to design 
molecular catalysts that employ earth-abundant materials, to match up to or surpass the catalytic activity of Pt. 
Because high electric conductivity and good catalytic activity are the two main requirements for an efficient 
electrocatalyst, screening the alternatives among the abundant conductive materials like conductive oxides is a 
reasonable way. However, according to our previous theoretical research 7 , some conductive oxides such as 
indium oxide (ln 2 0 3 ), stannic oxide (Sn0 2 ) and zinc oxide (ZnO) are unfortunately proven to be electrocata- 
lytically inactive, due to their low adsorption energies of reactant molecules and scarce active sites which are 
necessary for adsorption of the reactants, bond-breaking and bond-formation, and desorption of the products in 
the heterogeneous catalytic process 8 . To deliver these intrinsically inactive but abundant materials into active 
catalysts, one needs to control the atomic-scale surface structure to preferentially expose a greater fraction of the 
active sites 9 or exhibit a higher adsorption energy, which is highly challenging but desirable. 

Along the way, recent studies have shown that doping the carbon nanotubes and graphene with heteroatoms or 
defects sites can introduce electrocatalytic active sites for oxygen reduction in fuel cells 1,10,11 and for l 3 ~ reduction 
in DSCs 6,12 . Similarly, nitrogen-doped indium oxide (N-In 2 0 3 ) can be photocatalytically active for water split- 
ting 13,14 . More remarkably, the nitrate-type species in interstitial sites have been identified as an important role in 
the photoelectrochemical improvement of N-In 2 0 3 14 . However, such design strategies have rarely been based on 
molecular-level understanding of the reaction mechanism. In addition, how the heteroatoms transform the metal 
oxide with scarce active sites into high-efficiency catalysts with sufficient active sites has been unexplored. Within 
the density functional theory (DFT) framework, the adsorption energy of I atom (£ a d') on the common (110) 
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surface of ln 2 0 3 is found to be increased significantly when the 
interstitial N atoms are imported, which predicts that interstitial 
N-In 2 0 3 could be catalytically active. 

Herein, directed by this prediction, we describe our efforts to turn 
abundant ln 2 0 3 into a superior electrocatalyst by atomic-scale dop- 
ing. We accomplish this through a facile one-step solution strategy to 
introduce trace amount of interstitial N atoms to ln 2 0 3 nanocrystals 
via decompostition of indium nitrate, which serves as both reactant 
and dopant source. To the best of our knowledge, this is the first 
report of interstitial N-In 2 0 3 nanocrystals synthesized by a wet 
chemical technique. Further, the resultant N-In 2 0 3 nanocrystals 
supported by fluorine-doped tin oxide (FTO) glass substrates were 
applied as counter electrode (CE) in DSCs for the first time. 
Amazingly, the DSCs with N-In 2 0 3 CE achieved a power conversion 
efficiency (17) of 7.78% versus 7.17% for Pt under the same condi- 
tions, whereas the undoped ln 2 0 3 shows a 17 of only 0.96%. 
Furthermore, in order to understand how the trace nitrogen atoms 
turn the inactive ln 2 0 3 into a superior catalyst over Pt for I 3 ~ reduc- 
tion, the catalytic activity of interstitial N-In 2 0 3 towards I 3 ~ reduc- 
tion was investigated by means of DFT calculations in detail. 

Results 

To explore the effect of interstitial N-doping on the catalytic activity 
of ln 2 0 3 toward I 3 ~ reduction reaction, we firstly performed DFT 
calculation to investigate the doping configuration in ln 2 0 3 and its 
adsorption toward I atom, aiming at providing a preliminary under- 
standing for experimental synthesis. In pure ln 2 0 3 bulk, there exist 
two kinds of six-coordinated In atoms, forming a relatively regular 
and distorted octahedron structure, in which almost all the (at least 
seven kinds) possible initial structures for interstitial N-doping were 
optimized to identify the most stable one. The most stable config- 
uration achieved is shown in Figure la. One can see that the inserted 
N atom prefers to bond with two lattice O atoms and an In atom, 
forming a local N0 2 5 ~ configuration with a bond length of N-O of ~ 
1.36 A, in which the two oxygen atoms deviate evidently from their 
original lattice site in ln 2 0 3 and the number of O-In bond is reduced 
from 4 to 2. Bader charge analysis indicates that these two O atoms 

are charged with 0.77 e that is evidently less than the other O 

atoms (' — 1.26 e), and N exhibits an oxidized state (denoted as N 5+ 
with 5 = 0.14 e). Based on this stable configuration, ln 2 0 3 (l 10), one 
of the most exposed surfaces of ln 2 0 3 , was modelled with interstitial 
N dopant at the near-surface (see Figure lb and c). It shows that the 
doped N atom tends to be located at the subsurface and bind with the 
surface five-coordinated In atom at a distance of 2.30 A, thus it can 
be expected that the binding ability of this kind of surface In atom 
toward I or I 2 could be tuned due to the change of coordination 
environment. We calculated the adsorption of I atom at this five- 
coordinated surface In atom at the gas/solid interface, and the 
adsorption energy (relative to half an I 2 molecule) was calculated 
to be —0.54 eV for N doped In 2 O 3 (110), which is much enhanced 
compared to that (0.45 eV) of the pure In 2 O 3 (110) (see Figure lb). 



The adsorption of I atom at the interfaces of CH 3 CN/In 2 O 3 (110) 
with or without N doping was also calculated, which increase to 
— 0.94 and 0.16 eV, respectively, but give the same trend. Based on 
the thermodynamics model verified previously 7,15 , demonstrating 
that the £ a a : at the CH 3 CN/electrode interface can serve as a 
descriptor for the iodine reduction activity, and it should usually 
range from around —0.33 to —1.20 eV for good catalysts. Thermo- 
dynamically, when £ ad : is evidently weak, the dissociation of I 2 mole- 
cule would hardly proceed; while £ a< j : is too strong, the overall 
activity can be limited by I* removal into I~(sol) owing to the large 
energy requirement. It is clear that the adsorption energy of I atom 
on ln 2 0 3 resulting from the interstitial N doping is well within 
required thermodynamic range for a good catalyst to catalyze triio- 
dide reduction, suggesting the potentially enhanced catalytic activity 
of N-In 2 0 3 . 

Directed by the theoretical estimation, N-In 2 0 3 nanocrystals were 
synthesized by a smart and facile one-pot solution process (see the 
Methods). The scanning electron microscopic energy- dispersive x- 
ray (SEM-EDX) analysis in Figure SI (shown in SI) indicates that 
trace amount of N atoms of about 5.46 at.% were doped. In order to 
examine chemical changes caused by the doping with N atoms and 
the doping position of N atoms in ln 2 0 3 nanocrystals, the X-ray 
photoelectron spectroscopic (XPS) measurements were carried out. 
The XPS spectra of N-In 2 0 3 sample recorded is shown in Figure 2a-b 
for N 1 s and O 1 s core levels, and the survey spectrum and high- 
resolution In 3d XPS are shown in Figure S2-S3 in SI. Nitrogen 1 s 
core level shows a single peak at 399.3 eV, which appears at higher 
binding energy (BEs) than the typical energy for InN (397 eV 16 , or 
396.5 eV 17 ) and at lower BEs than free or adsorbed N species on 
metal oxides. According to the literature about N-doped 
ln 2 0 3 13,14,18 or N-doped Ti0 2 19 , nitrogen atoms substitutionally 
doped into the ln 2 0 3 lattice (In-N-In linkages) own the characteristic 
peak at approximate 397 eV, which does not appear in the XPS of our 
sample, and on the contrary, the peak at 399.3 eV in our sample can 
be attributed to a characteristic peak of interstitial N, being NO x 
located in an interstitial position. The N-In 2 0 3 sample shows a single 
peak at 444.0 eV for In 3d 5/2, which is attributed to ln 2 0 3 according 
to the NIST XPS database (444.30 to 446.7 eV). This small lower BE 
shift may be due to a interstitial doping 13 . Gaussian fitting of the O 1 s 
region shows two peaks at 529.4 and 531.0 eV. The former peak can 
be attributed to the In-O bond in ln 2 0 3 13,18 , and the peak at 531.0 eV 
is the characteristic of a non-stoichiometric oxide, oxynitride, or 
adsorbed water 17,18 . It is worth noting that the peak at 531.0 eV over- 
whelms the peak at 529.4 eV of pristine ln 2 0 3 . This result indicate 
that the interstitial N enables the sample adsorb the foreign hydro- 
xyls more easily 18 , which is consistent with our theory prediction of 
the increasing at the acetonitrile/electrode interface. The mor- 
phology and crystallinity of as-prepared N-In 2 0 3 nanocrystals were 
analyzed using transmission electron microscopy (TEM) and high- 
resolution TEM (HRTEM). Figure 2c shows the TEM image of quasi- 
monodisperse N-In 2 0 3 nanocrystals and the nanocrystal diameter 




Figure 1 | Structures of interstitial N doping and I atom adsorption, (a) The most stable bulk structure of interstitial N-doped ln 2 0 3 (N-In 2 0 3 ); (b) top 
view of the optimized surface structure of N-In 2 0 3 (110); (c) side view of the optimized I adsorption structure on N-In 2 0 3 (110). Red, brown, blue and 
purple balls represent O, In, N and I, respectively. Marked in circle shows the local configuration of N0 2 6 ~ resulting from the interstitial N doping, and 
dark brown balls indicate the distorted lattice O accompanying the N doping. 
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Figure 2 | Atoms distribution and morphological analysis of N-In 2 C>3 nanocrystals. (a) high-resolution N 1 s XPS; (b) high-resolution O 1 s XPS; The 
open circles are the raw data of the X-ray photoelectron spectra, and the blue and red/purple/green lines represent the base line and fitted lines, 
respectively, (c) representative morphologies and structures of the N-In 2 03 nanocrystals in TEM; (d) HRTEM images. 



was approximately (9 ± 2) nm. From Figure 2d, it can be seen that 
the N-In 2 0 3 nanocrystal is mainly bounded by {440} facets. As 
shown in Figure S4 in SI, all X-ray diffraction (XRD) peaks of the 
synthesized N-In 2 0 3 can be indexed to the lattice planes of body- 
centered cubic ln 2 0 3 (JPCDs no. 71-2194), which is consistent with 
the TEM result. It is worth noting that compared to conventional 
costly preparation of irregular N-In 2 0 3 powders via high-temper- 
ature calcinating, the way to produce N-In 2 0 3 nanocrystals with 
small grain size and faceted morphology via a solution-base 
approach would lower the manufacturing cost and increase the cata- 
lytic activity of the material. 

The performances of N-In 2 0 3 on I 3 ~ reduction and further on 
photovoltaic conversion efficiency of DSCs were characterized in the 
following parts. For comparison, platinum and undoped ln 2 0 3 
nanooctahedron synthesized according to the literature (see SI in 
details, SEM shown in Figure S5) 2 ", were also used as CEs in DSCs 
constructed under the same condition in control experiments. 
Figure 3a shows the photocurrent-voltage (/-V) curves of three 
DSCs fabricated with N-In 2 0 3 , Pt and undoped ln 2 0 3 nanooctahe- 
dron electrodes as CEs. The detailed photovoltaic parameters from 
the /-V curves are summarized in Table 1. As shown in Table 1, the 
device equipped with N-In 2 0 3 yields a remarkable ;/ of 7.78%, which 
is superior to that of the device with a Pt CE, far better than that of 
undoped ln 2 0 3 nanooctahedron and / sc increased in the same order. 
These interesting results demonstrate that the interstitial N atoms 



surely turn ln 2 0 3 into a superior catalyst surpassing Pt for I 3 ~ reduc- 
tion, which is in accordance with our prediction results. To better 
understand the improved DSC performance for N-In 2 0 3 CE, elec- 
trochemical impedance spectra (EIS) experiments of the symmetrical 
cells fabricated with two identical N-In 2 0 3 and Pt CEs were carried 
out at 20°C, and the results are shown in Figure 3b, while the EIS of 
undoped ln 2 0 3 nanooctahedron was given in SI (Figure S6). The EIS 
parameters were determined by fitting the impedance spectra using 
the ZSimpWin software via the circuit model inserted in Figure 3b, 
and were summarized in Table 1. The high-frequency intercept on 
the real axis represents the series resistance (R s ), which is mainly 
composed of the bulk resistance of CEs materials, resistance of 
FTO glass substrate, and contact resistance, etc. It can be seen that 
R s of N-In 2 0 3 (7.03 fi) is smaller than that of Pt (9.64 Q), which 
should be benefit to the increase of the / sc of N-In 2 0 3 . While the 
lower R s of N-In 2 0 3 than that of undoped ln 2 0 3 (9.52 Q) indicates 
that the introduction of the interstitial N atoms increases the electric 
conductivity of ln 2 0 3 . The low-frequency region (the right semi- 
circle) reflects the diffusion in the electrolyte, which is not closely 
related to catalysis and hence is not the focus of attention in this 
paper. The semicircle in the high frequency region (the left one) 
arises from the parallel connection of coupling between the constant 
phase angle element (CPE) at the CE/electrolyte interface and the 
charge-transfer resistance (.R ct ) of CE, which is directly proportional 
to the diameter of the left semicircle. R a is a measure of the ease of 
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Figure 3 | DSCs performances and electrochemical analysis. 

(a) Photocurrent- voltage (I-V) curves of DSCs fabricated with Pt, N-In 2 0 3 
and ln 2 0 3 as counter electrodes and (b) electrochemical impedance 
spectra of the symmetrical cells fabricated with two identical counter 
electrodes, and the insert gives the equivalent circuit where W( so i) 
represents Warburg parameter. 

electron exchange between the electrode and the electrolyte and thus 
varies inversely with the I 3 ~ reduction activity of the CEs 21 . As shown 
in Table 1, the R ct of N-In 2 0 3 is 5.51 Q, significantly smaller than 
that of Pt (8.65 Q) and much smaller than that of undoped ln 2 0 3 
nanooctahedron (379.4 Q), which means that the N-In 2 0 3 has a 
higher catalytic activity for I 3 ~ reduction than that of Pt. While the 
so high _R ct of undoped ln 2 0 3 nanooctahedron indicates that it nearly 
has no catalytic activity for I 3 ~ reduction. It also can be seen that the 



CPE of N-In 2 0 3 is lower than that of Pt, which means that N-In 2 0 3 
has a smaller active surface area 22 24 . But, it can be obtained that the 
low R s and high catalytic activity (_R ct ) of N-In 2 0 3 make up its small 
active surface area, resulting in the high / sc and r\ of DSCs with N- 
ln 2 0 3 . In order to further elucidate the catalytic activities of different 
CEs, we also performed the cyclic voltammograms (CVs) of the 
ir/l - redox couple on both the N-In 2 0 3 and Pt electrodes. For the 
both electrodes, two pairs of oxidation and reduction peaks are well 
resolved, as shown in Figure S7. The left and right pairs are described 
by equations (1) and (2), respectively. 

I--f-2e~<-»3I~ (1) 



3I 2 + 2e-^2FT (2) 

Since the CE of a DSC is responsible for I 3 ~ reduction, we only focus 
on the left cathodic peak corresponding to I 3 ~ reduction. The high 
peak current and large peak area of the left cathodic peak reveal that 
N-In 2 0 3 is a significant efficient electrochemical catalyst for I 3 ~ 
reduction, even better than Pt electrode 6,25,26 . To further examine 
the interfacial charge-transfer properties of the I 3 ~/I~ couple on 
the electrode surface, Tafel polarization measurements were carried 
out with dummy cells, similar to the one used in EIS measurements, 
and the results are shown in Figure S8 in SI. From the Tafel zone at 
middle overpotential (with a sharp slope), we can obtain the 
exchange current density (/ 0 ), which can also be obtained from EIS 
parameters via equation (3). 



where R ct is determined by EIS data, and R, T, n and F have their usual 
significance, which means that high J 0 means low _R ct and high cata- 
lytic activity. The anodic and cathodic branches of the J-U plot 
(Figure S8 in SI) show a largest slop for the N-In 2 0 3 , indicating a 
highest / 0 on this electrode in terms of the Tafel equation 27,28 , which 
means that N-In 2 0 3 has a higher catalytic activity. These experi- 
mental results consistently show that the interstitial N atoms surely 
turn the intrinsically inactive ln 2 0 3 into an efficient electrocatalyst. 

Further, to elucidate the origin of high catalytic activity of N-In 2 0 3 
and understand the role of the interstitial N in the performance of 
catalyst, a molecular level investigation into I 3 ~ reduction reaction at 
the CH 3 CN/N-In 2 O 3 (110) interface was conducted by virtue of DFT 
calculations. Herein, a full reaction energetic involving thermodyn- 
amics and kinetic are considered to demonstrate the catalytic activity 
of N-In 2 0 3 toward I 3 ~ reduction. With respect to the I 3 ~ reduction 
reaction (I 3 ~( so i) +2e"^ 3I~( so i)) occurring on the CE, the general 
consensus of the reaction mechanism can be described in the refer- 
ence 7,29 (see details in SI), which shows that the I 2 reduction reaction 
(IRR) occurring at the liquid-solid interface, can determine the over- 
all electrocatalytic activity owing the fast solution reaction (I 3 ~( so i) <-» 
l2(soi) + I (sol))- Hence we focused our studies on the two elementary 
reaction steps involved in I 2 reduction reaction at the CH 3 CN/elec- 
trode interface, containing I 2 ( so i) + 2* — > 21*; I* + e~ — » T~( so i)> to 
explore the activity of ln 2 0 3 with or without N doping. It is found 
that on N-In 2 O 3 (110) surface I 2 molecule can readily dissociate into 
two I* atoms without an obvious dissociation barrier (-E a dls )> resulting 
in I* sitting preferably on the top of In atom at a distance (d(i n -i)) of ~ 
2.93 A, while it is hard on pure In 2 O 3 (110). The adsorbed iodine 
atom (I*) can receive one electron and then desorb into solution in 
form of I with a much elongated d( Tn _i) of 4.20 A at the transition 
state (TS) structure, being surrounded by CH 3 CN molecules with H 
atoms in methyl pointing to it. The charge density differences show 
(Figure 4c) that the electrons accumulate at I atom in TS (Bader: 
0.69 e), resulting from the electron depletion at the surface In atom 
and the surface adsorbed CH 3 CN molecules, and there is no orbital 



Table 1 | performance of the DSCs using Pt and N-ln 2 03 as CEs [al 
and EIS parameters of the symmetrical cells fabricated with two 
identical counter electrodes [bl 

Samples i/mA cm" 2 V oc /mV FF RjQ RjQ. CPE/uF 

Pt 18.07 749 0.54 7.17 9.64 8.65 7.66 

N-ln 2 0 3 19.33 750 0.54 7.78 7.03 5.51 5.41 
ln 2 0 3 11.13 336 0.26 0.96 9.52 379.4 4.56 

short-circuit current density; 
VqcI open-circuit voltage; FF: fill factor; //: energy conversion efficiency. 
[b| For spectra, see Figure 3(b). 

R s : series resistance; R^: charge-transfer resistance at the CE/electrolyte interface; CPE: constant 
phase-angle element at the electrolyte/CE interface. 
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overlap between thep-orbital of ion-like I andp-orbital of surface In 
atom, indicating a typical ionic bond of In-I. Figure 4a shows the 
reaction standard Gibbs free energy profile of IRR at the CH 3 CN/N- 
In 2 O 3 (110) interface, in which the energy profiles on two typical 
surfaces of Pt electrode, i.e. flat Pt(lll) and stepped Pt(411), were 
also given as the references 7,25 . From the profile of N-In 2 0 3 , one can 
see that the Gibbs free energies of all elemental reactions at CH 3 CN/ 
N-In 2 0 3 (l 10) interface decrease step by step from the initial state to 
final state, indicating a considerable catalytic activity. Quantitatively, 
the adsorption energy of I atom at CH 3 CN/N-In 2 0 3 (l 10) interface is 
between Pt(411) and Pt(lll), giving a —0.94 eV, while the corres- 
ponding desorption barrier of I* is only 0.28 eV, lower than those at 
the CH 3 CN/Pt(lll) and CH 3 CN/Pt(411) interfaces (0.39 eV and 
0.63 eV, respectively). Considering that I 2 molecule can readily dis- 
sociate on all these three surfaces, the dissociation of I 2 molecule 
could be expected to be fast and may be in equilibrium approxi- 
mately, and one can therefore speculate that N-In 2 0 3 can exhibit a 
good catalytic activity toward IRR owing to the lower £ a des , even 
better than Pt(lll). Furthermore, we estimated the equilibrium 
exchange current (i 0 ) of IRR using the micro-kinetic analysis to 
compare their activity trend. For these three surfaces, the i 0 of the 
IRR can be determined from equation ( 1 0) in SI, which quantitatively 
indicates the key role of the desorption barrier of I* (E a des ); the 
smaller the £ a d<:s , the higher the catalytic activity. Accordingly, it 
can be concluded that their catalytic activities toward IRR follow in 
the order of N-In 2 O 3 (110) > Pt(lll) > Pt(411), being consistent 
with the experimental R ct values, CV and Tafel measurements. 

At this point, one question arises: why does interstitial N doping 
give rise to the enhancement of I* adsorption on ln 2 0 3 (l 10) surface? 
Geometrically, the surface five-coordinated In atom on In 2 O 3 (110) 
surface strongly bonds with four surface lattice O (O s ) and one sub- 
surface O (O ss ). After interstitial N doping at the subsurface, this 
surface In atom would form a relatively weak In-N bond compared 
to In-0 bond, reflected in a longer In-N bond length (2.37 A vs. 
2.16 A for In-O ss bond). On account of the bond conversation, the 
enhanced binding ability of the surface In atom resulting from N 
doping is reasonable. Electronically, as shown in Figure 4b, the I-In 
bond at the I* adsorption state exhibits evident characteristic of 
ionic bond, and thus Bader charge analysis is conducted to under- 
stand the N doping effect (see details in Table S2 in SI). It shows 
that the electrons of the ln 2 0 3 system owing to N doping would be 



redistributed, probably rationalized by the obvious structural distor- 
tion in forming a local N0 2 6 ~ doping configuration. Compared with 
the pure ln 2 0 3 , the total electron of all the O atoms would decrease by 
1.04 e, and the two O atoms binding directly with N atom decrease 
most, each giving ~ 0.49 e; while the total electron number of all the 
In atoms increases by 1.14 e, and the increased electrons were delo- 
calized at all of In atoms with a little more distributing at the surface 
In atoms. Upon adsorption on N-In 2 0 3 (l 10), I atom can gain 0.47 e, 
which mainly results from the In atoms (0.28 e) and the interstitial N 
dopant (0. 10 e), giving rising to a relatively strong adsorption energy 
toward I atom ( — 0.54 eV). By contrast, as In 3+ in pure ln 2 0 3 hardly 
further provide electron and there are no extra electrons on them, I* 
can only gain few electrons (0.29 e) from the O atoms (0.28 e) upon 
adsorption on pure In 2 O 3 (110) surface, resulting in a weak ionic 
bond and low adsorption energy (0.45 eV). In general, the enhance- 
ment of adsorption ability of ln 2 0 3 due to interstitial N doping could 
mainly be ascribed to the electron redistribution of ln 2 0 3 , in which 
more extra electrons can distribute on In atom to facilitate I* adsorp- 
tion, as well as the N atom itself would provide a small number of 
electrons in forming a local N0 2 6 ~ doping configuration. 

As a whole, both experimental measurements and first-principles 
calculations demonstrate that the performance of ln 2 0 3 catalyst in 
catalyzing I 3 ~ reduction can be optimized by doping N atom, and the 
catalytic activity of N-In 2 0 3 can be even better than Pt(l 1 1) owing to 
its increased I* adsorption energy and lower desorption barrier. 

Discussion 

In summary, for the first time, we have prepared interstitial nitrogen- 
doped ln 2 0 3 nanocrystals via a low cost and facile one- step solution 
strategy, and demonstrated that the small amounts of interstitial 
nitrogen can transform the intrinsically inactive ln 2 0 3 into a high- 
efficiency catalyst in DSCs, leading to a high /; of 7.78%. Beyond this 
application, this work demonstrates that a cheap and abundant semi- 
conductor with high electric conductivity such as ln 2 0 3 can be 
turned into an active catalyst by controlled interstitial doping of N 
atoms via a novel nanotechnology of a simple solution-based strat- 
egy, and opens a new pathway to prepare effective Pt-free electro- 
catalysts. More importantly, the exciting discovery of turning the 
commonly used TCO in DSCs into a difunctional electrode material 
with both transparent conductive function and catalytic activity 
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Figure 4 | Electronic structures of intermediates and energy profiles, (a) The standard Gibbs free energy profile of IRR at various systems of Pt( 1 1 1 ), 
Pt(411) and N-In 2 O 3 (110); (b) and (c) indicate the structure of I adsorption and desorption transition state at the CH 3 CN/N-In 2 O 3 (110) interface, 
respectively, in which 3D isosurfaces of charge density difference before and after I atomic interacting with the electrode are also illustrated 
(yellow indicates the electronic accumulation and light blue for electronic depletion). 
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means that except for decreasing the cost, the device structure and 
processing techniques of DSCs can be simplified in future. 

Methods 

Synthesis of N-In 2 0 3 nanocrystals. In a typical synthesis of interstitial N-In 2 0 3 
nanocrystals, 0.18 g In(N0 3 ) 3 - 4.5H 2 0 was dissolved in the mixture of 5 mL 1- 
octadecylamine and 5 mL 1-octadecene (at 65' C). The solution was degaseed for at 
least 2 hours and heated at 150 "C for 1 h under magnetic stirring. Then the solution 
temperature was raised to 280°C in 10 min and kept for 30 min. The product was 
collected by centrifugation, washed with chloroform and alcohol several times. 

Materials characterizations. The morphology and structure of the samples were 
characterized by high- resolution transmission electron microscopy (TEM; Tecnai 30, 
300 kV). XRD patterns were recorded on the dry membrane of aqueous dispersions of 
the N-In 2 0 3 nanocrystals on a glass wafer by a Bruker D8 Advanced Diffractometer 
(Japan) using Cu K a radiation (X = 1.5406 A) at 40 kV. Binding elements were 
analyzed with X-ray photoelectron spectroscopy (XPS, Kratos Axis Ultra DLD) and 
CASAXPS software was used to analyze the XPS data. All binding energies were 
referenced to the Cls peak (284.8 eV) arising from surface hydrocarbons (or possible 
adventitious hydrocarbon). The photocurrent-voltage performance of the DSCs was 
measured with a Keithley digital source meter (Keithley 2400, USA) and simulated 
under AM 1.5 illumination. The power of the simulated light was calibrated to 
100 mW'cm" 2 using a Newport Oriel PV reference cell system (model 91150 V). All 
the electrochemical experiments were conducted by using a computer- controlled 
electrochemical workstation (Parstat 2273, Princeton) in the dark. The measured 
frequency of EIS experiments ranged from 100 mHz to 1 MHz with an AC 
modulation signal of 10 mV and bias DC voltage of 0.60 V. The spectra were fitted by 
the ZSimpWin software. Tafel polarization measurements were carried out in a 
dummy cell with a scanning rate at 50 mV s~\ Cyclic voltammetry (CV) was 
conducted in a three- electrode system in an acetonitrile solution of 0.1 M LiC10 4 , 
10 raM LSI, and 1 mM I 2 at a scanning rate of 20 mV s" 1 . The working electrodes were 
N-In 2 0 3 nanocrystals and Pt loaded FTO, exposing an area of about 0.25 cm 2 . The Ag/ 
Ag + couple was employed as reference electrode, and the CE was a platinum foil. 

Theoretical calculation. All the spin-polarized calculations were performed with 
Perdew-Burke-Ernzerhof (PBE) functional 30 within the generalized gradient 
approximation using the VASP code 31,32 . The project-augmented wave (PAW) 
method was used to represent the core-valence electron interaction 33,34 . To model the 
In 2 O 3 (110) surface, a seven-layer p(l X 1) slab (14.406 X 10.187 A 2 ) 35 " 40 
corresponding to 28 ln 2 0 3 units cell ( 1 40 atoms) was used, in which a vacuum layer of 
15 A was applied. Because of the large size of the supercell, k-point sampling was 
restricted to the V point only. For total energy calculations, the valence electronic 
states were expanded in plane wave basis sets with an energy cutoff of 450 eV, and the 
force convergence criterion in structural is set to be 0.05 eV/A. The transition states 
were searched using a constrained optimization scheme 41 , and were verified when (i) 
all forces on atoms vanish; and (ii) the total energy is a maximum along the reaction 
coordination but a minimum with respect to the rest of the degrees of freedom. The 
force threshold for the transition state search was 0.05 eV/A. To simulate the 
interface solution environment, several acetonitrile layers were introduced in the 
surface slab with a density of 0.79 g cm -3 . The adsorption energy of I (E a d) is defined 
as: £ a d J = £( I/interface) — ^(interface) + 1/2£(I 2 ), where £( interface), E(l 2 ) and £(1/ 
interface) are the energies of the liquid/electrode interface, I 2 in the gas phase and I 
adsorbed on the liquid/electrode interface, respectively. In addition, please see details 
for the relevant free energy calculation of each elementary reaction step. 
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